SFMAR 


Systems  Center 
PACIFIC 


TECHNICAL  REPORT  3056 
February  2017 


Quasi-static  Design  of  Electrically  Small 

Ultra-Wideband  Antennas 


Thomas  O.  Jones  III 


Approved  for  public  release. 


SSC  Pacific 
San  Diego,  CA  92152-5001 


SSC  Pacific 

San  Diego,  California  92152-5001 


G.  M.  Bonitz,  CAPT,  USN  C.  A.  Keeney 

Commanding  Officer  Executive  Director 


ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  performed  by  the  Electromagnetics  Technology  Branch 
(Code  52260)  of  the  Systems  of  Systems  (SoS)  and  Platform  Design  Division  (Code  52200),  Space 
and  Naval  Warfare  Systems  Center  Pacific  (SSC  Pacific),  San  Diego,  CA. 


Released  by 
M.  Osburn,  Head 

Electromagnetics  Technology  Branch 


Under  authority  of 
J.  McGee,  Head 
Systems  of  Systems  (SoS)  and 
Platform  Design  Division 


This  is  a  work  of  the  United  States  Government  and  therefore  is  not  copyrighted.  This  work  may 
be  copied  and  disseminated  without  restriction. 


The  citation  of  trade  names  and  names  of  names  of  manufacturers  is  not  to  be  construed  as  official 
government  endorsement  or  approval  of  commercial  products  or  services  referenced  in  this  report. 

CST  Microwave  Studio®  is  a  registered  trademark  of  Dassault  Systems. 


SB 


EXECUTIVE  SUMMARY 


Ultra  wideband  antennas  (UWB)  can  be  designed  with  electrostatic  methods  using  a  generalization 
of  the  [1]  asymptotic  conical  dipole  (ACD).  The  electrostatic  methods  allow  the  antenna  (J-factor  to 
be  calculated  for  a  large  number  of  antenna  designs  without  numerically  solving  Maxwell's 
Equations.  The  ACD  uses  a  constant  line  charge  distribution  and  image  line  charge  distribution  (both 
on  the  3-axis)  to  generate  equipotential  surfaces.  Each  equipotential  surface  represents  an  ACD 
antenna  design  with  a  different  height.  In  the  Quasi-static  Antenna  Design  Algorithm  [2,  3,  4,  5,  6] 
the  antenna  designs  fit  within  a  1-meter  sphere.  Both  the  charge  distribution  and  charge  distribution 
geometry  are  varied  to  generate  a  range  of  antenna  shapes.  The  Q-factor  ratio,  radiation  resistance, 
and  capacitance  design  curves  are  calculated  from  the  geometry  and  charge  distribution.  Each  charge 
distribution  and  geometry  represents  a  unique  antenna  shape  and  surface  charge  density.  Solutions 
are  selected  from  the  design  curves,  then  the  antenna  shape  is  calculated;  the  last  step  is  numerical 
electromagnetic  modeling  of  the  antenna.  For  the  UWB  antenna  designs,  the  charge  distribution  is  on 
a  cone  and  an  image  cone.  The  charge  distribution  on  the  cone  and  the  image  cone  is  selected  to  give 
a  conical  feed  point.  The  cone  half  angle  and  cone  length  are  the  geometric  parameters  used  to 
generate  the  Q-factor  ratio,  radiation  resistance,  and  capacitance  design  curves.  The  minimum 
Q-factor,  radiation  resistance,  and  capacitance  for  the  UWB  antenna  designs  are  plotted  as  a  function 
of  the  charge  distribution  half  cone  angle.  In  this  plot,  the  radiation  resistance  and  capacitance  are 
discontinuous.  This  discontinuity  is  caused  by  the  absolute  and  local  ratio  minimum  in  the  Q-factor 
ratio  design  curve.  The  discontinuity  occurs  when  two  different  charge  distribution  cone  geometries 
give  the  same  Q-factor  ratio.  Three  antennas,  50-,  30-,  and  (the  minimum  Q-factor  ratio)  32-Q.  input 
impedance  designs  that  are  numerically  modeled  from  30  MHz  to  2.5  GHz.  The  antenna  impedance 
and  voltage  standing  wave  ratio  (VSWR)  is  calculated.  In  all  models,  the  VSWR  drops  sharply  to 
about  1.2:1,  then  drops  to  1.1:1  for  nearly  all  of  the  frequency  range.  The  Smith  chart,  time  domain 
feed  point  input,  and  feed  point  reflection  are  plotted.  The  minimum  (^-factor  ratio  design  has  a 
3:1  VSWR  at  42.67  MHz  or  an  antenna  height  of  one-seventh  of  a  wavelength  at  the  lowest 
operational  frequency  and  8.34  wavelengths  at  the  highest  operational  frequency. 
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1.  INTRODUCTION 


The  ultra  wideband  (UWB)  antenna  shape  determines  both  the  lowest  operating  frequency  and 
(^-factor  ratio.  For  electrically  small  antennas,  the  quasi-static  antenna  design  algorithm  can  be  used 
to  design  very  low  Q- factor  ratio  antennas  with  high  bandwidth  [4,  5].  The  (^-factor  ratio  is  the 
coefficient  of  the  largest  term  in  Chu's  limit,  the  right  side  of  Equation  (1)  [7].  Lowering  the  Q-factor 
ratio  reduces  the  (^-factor  for  the  electrically  small  part  of  the  frequency  range.  The  UWB  antenna 
shape  in  the  quasi-static  antenna  design  algorithm  is  a  cone  with  an  attached  bulb.  The  lowest 
operational  frequency  is  reduced  without  increasing  the  voltage  standing  wave  ratio  (VSWR)  at 
higher  frequencies.  This  is  a  physics-based  algorithmic  approach  that  avoids  the  typical  cut- and- try 
method. 

In  1948,  Chu  [7]  established  a  lower  bound  for  the  (^-factor  (inverse  bandwidth)  of  an  antenna. 
Chu's  limit  is  calculated  from  the  energy  stored  in  the  electromagnetic  field.  The  antenna  is 
enclosed  within  a  sphere  with  a  radius  a;  only  the  stored  energy  outside  the  sphere  is  included  in 
the  limit.  The  energy  inside  the  sphere  is  zero.  For  top-loaded  monopoles,  the  energy  is 
electrostatic  below  the  antenna's  resonant  frequency  and  magnetic  above  the  antenna's  resonant 
frequency.  At  resonance,  the  electric  and  magnetic  energy  are  equal  but  90°  out  of  phase.  The  top- 
loaded  monopoles  have  a  large  reactance  and  stored  magnetic  energy  above  resonance.  On  the 
other  hand,  the  UWB  antennas  have  small  reactance  and  stored  magnetic  energy  above  resonance. 

In  UWB  antennas,  the  (^-factor  is  —  1  above  the  lowest  operational  frequency. 

An  infinite  bi-conical  (transmission  line)  is  a  UWB  antenna  with  a  constant  impedance;  the  cone 
angle  determines  the  impedance.  A  finite  bi-conical  antenna  was  patented  by  Oliver  Lodge  [8,  9]. 

In  his  design,  the  end  of  the  cone  reflects  the  outgoing  wave  back  into  the  feed  point.  Thus,  the 
antenna  is  no  longer  UWB  or  constant  impedance.  Reducing  this  reflection  from  the  top  of  the  cone 
is  a  longstanding  problem.  The  Schelkunoff  and  Friis  [10]  teardrop  antenna  (bi-cone  at  feed  point) 
improves  the  performance.  Kraus's  volcano  and  smoke  antenna  design  was  an  effort  to  eliminate  the 
reflection  from  the  antenna.  Paulson,  West,  Perger,  and  Kraus  [11]  state  that  the  volcano  and  smoke 
monopole  is  typically  one-fourth  a  wavelength  high.  Maloney  and  Smith  [12]  used  resistive  material 
in  the  cone  to  absorb  the  reflection  from  the  top  of  the  cone;  this  reduces  the  antenna  efficiency.  The 
Ice  Cream  Cone  (ICC)  monopole  designed  by  Simpson,  Pavlovic,  and  Olcan  [13]  does  not  use 
resistive  material.  The  combined  spherical  bulb  and  cone  has  no  reflection  from  the  feed  point.  The 
ICC  antenna  height  is  one-fifth  of  the  longest  operational  wavelength1.  For  a  1 -meter  monopole,  the 
lowest  operational  frequency  is  60  MHz.  Another  approach  to  the  UWB  antenna  design  is  the 
exponential  taper  of  Hale  Simonds  [14]  and  hyperbolic  cosine  taper  of  Thomas  O.  Jones  III  [15];  this 
design  reduces  the  width  of  the  antenna,  which  implies  that  the  bulb  shape  can  be  non-spherical  at 
high  frequencies. 

The  stored  energy  in  an  UWB  antenna  is  analyzed  differently.  The  cone  is  a  transmission  line  that 
supports  an  outgoing,  freely  propagating  wave;  the  outgoing  wave  does  not  represent  stored  energy. 
The  transmission  line  is  terminated  with  a  bulb.  The  bulb  reflects  some  of  the  radio-frequency  (RF) 
energy  back  into  the  transmission  line  to  the  feed  point.  This  reflected  energy,  an  incoming  wave, 
is  included  in  the  stored  energy  and  in  the  Q-factor.  The  bulb  also  contributes  directly  to  the  stored 
energy.  The  bulb  shape  is  the  critical  factor  in  determining  the  Q-factor.  For  a  fixed  cone  angle,  the 
cone  plays  no  role  in  the  (^-factor.  This  report  presents  a  scientific  method  for  designing  low 
Q-factor  UWB  bulb  shapes. 


1  Personal  communication  with  Ted  Simpson. 
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The  Quasi-static  Antenna  Design  Algorithm  is  used  to  design  a  minimum  (J-factor  bulb  shape  with 
a  conical  feed  point.  The  algorithm  uses  electrostatic  methods  [2,  3,  4,  5,  6],  In  electrically  small 
antennas,  the  magnetic  fields  are  small  and  can  be  neglected,  which  permits  the  use  of  the  quasi¬ 
static  approximation  used  in  the  algorithm.  A  static  charge  distribution  is  used  to  generate 
equipotential  surfaces.  The  equipotential  surfaces  are  perfect  conductors  with  E\  =  0.  Each  equipo- 

tential  surface  represents  a  unique  antenna  surface  and  surface  charge  density.  For  an  antenna,  the 
charge  is  only  on  the  surface  of  the  conductor.  The  Quasi-static  Antenna  Design  Algorithm  requires 
the  charge  distribution  to  be  enclosed  by  the  equipotential  surfaces. 

The  antenna  (Q-factor  is  calculated  from  radiation  resistance,  /?  ,  the  capacitance  C,  the  angular 

frequency  co  =  2nf: 


Q 


l 


iCR 


Rad 


Q  Factor  ratio 


( ka )3 


(1) 


where  /  is  the  frequency.  The  (^-factor  ratio  (Q-Factor  ratio)  is  calculated  from  Q and  ka,  where  k  = 

2n/X,  where  X  is  the  wavelength,  and  a  is  the  enclosing  spheres  radius.  The  radiation  resistance  is 
calculated  from  the  effective  height  and  the  effective  height  is  calculated  from  charge  distribution. 
The  capacitance  is  calculated  from  the  peak  scalar  potentials  on  the  enclosing  sphere  [2,  3],  For  a 
fixed  DC  voltage  V,  higher  capacitance  increases  the  stored  energy  E  =  CV2/2.  For  a  fixed  charge  Q, 
increasing  the  capacitance  decreases  the  stored  energy  E  =  Q2/2C.  The  Quasi-static  Antenna  Design 
Algorithm  used  a  fixed  charge  on  the  monopole  arm.  The  Quasi-static  Antenna  Design  Algorithm 
was  used  to  compute  the  thick-disk-cap  monopole  [13,  14],  The  charge  disk  height  and  radius 
describe  the  geometry  of  the  thick-disk  charge  distribution.  The  charge  distribution  on  the  disk  uses 
five  multipoles.  The  thick-disk-cap  monopole  required  only  seven  unknowns  [4,  5,  6]  to  converge  to 
the  solution.  Above  the  thick-disk-cap  monopole  resonance,  this  design  has  a  37%  lower  Q-factor 
than  the  spherical-cap  monopole  design  (the  best  known  top-loaded  monopole  design).  The 
approximations  in  the  Quasi-static  Antenna  Design  Algorithm  is  discussed  by  Thomas  O.  Jones  III  in 
[5,  6], 

The  Quasi-static  Antenna  Design  Algorithm  requires  a  charge  distribution  that  generates  a  wide 
range  of  cone  and  bulb  shapes.  The  z-axis  line  charge  distribution  used  in  [2,  3]  does  not  model  the 
general  bulb  shape.  A  rotationally  symmetric  three-dimensional  conical  charge  distribution  models 
the  bulb  shape  and  the  conical  feed  point.  Figure  1  is  a  MATFAB®  plot  of  the  50-  Cl  input 
impedance  UWB  antenna  with  a  120°  cutaway  showing  the  conical  charge  distribution.  The  color  is 
the  magnitude  of  the  surface  charge  density/m2  on  the  antenna  and  cone.  The  surface  charge 
density/m2  is  smaller  on  the  antenna  than  the  cone.  The  charge  on  the  antenna  is  spread  over  a  larger 
surface  area;  this  significantly  reduces  the  surface  charge  density.  In  addition,  the  image  cone 
reduces  the  electric  field  and  surface  charge  density  above  the  bulb.  The  image  cone  increases  the 
electric  field  and  surface  charge  density  near  the  feed  point  of  the  antenna. 

Section  2  calculates  the  effective  height  for  the  conical  charge  distribution.  The  charge  distribution 
is  selected  to  give  a  very  good  approximation  of  a  50-Q  cone.  This  conical  charge  distribution 
reduces  to  a  line  charge  on  the  z-axis  (cone  half  angle  ip  =  0).  The  normalized  radiation  resistance, 
RRJ(ka)2  is  computed  as  a  function  of  charge-distribution-cone  half  angle  and  fractional  cone  length 
(length  /a).  Section  3  develops  a  numerical  method  to  calculate  the  scalar  potential  for  the  general 
conical  charge  distribution.  The  numerical  method  uses  2N  wires  to  model  the  charge  distribution 
[14,  pp.  43-56  for  the  cylinder].  The  error  in  this  numerical  model  is  a  very  small,  high-order 
multipole  moment  in  the  scalar  potential.  The  high-order  multipole  moment  error  decreases  with 
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increasing  N.  The  capacitance  is  computed  as  a  function  of  charge-distribution-cone  half  angle  and 
fractional  cone  length. 


Figure  1 .  The  50-Q  input  impedance  antenna  shape  (equal-potential  surface)  with 
a  120°  slice  removed  to  show  the  conical  charge  distribution.  The  color  is  the 
charge/m2.  The  antenna  surface  charge  is  spread  over  a  much  larger  surface  area, 
thus  it  is  lower  than  the  surface  charge  on  the  cone. 

Section  4  calculates  the  (^-factor  ratio  as  a  function  of  the  charge-distribution-cone  half  angle  and 
fractional  cone  length.  The  (^-factor  ratio  design  curve  is  discussed  in  great  detail  in  Section  4.  For 
low-charge  distribution-cone  half  angles,  the  radiation  resistance  has  a  moderate  value.  At  a  critical 
charge  distribution-cone  half  angle,  the  radiation  resistance  jumps  to  a  much  higher  value.  At  higher 
charge  distribution-cone  half  angles,  the  radiation  resistance  drops,  the  capacitance  increases,  and  the 
Q-factor  ratio  drops  to  a  minimum  of  3.265.  The  antenna  parameters  are  calculated  for  several 
antenna  input  impedances:  50  £2,  30  Q,  and  for  the  minimum  Q-factor  antenna,  32  Q.  The  feed  point 
is  a  very  good  approximation  of  the  50-Q  input  impedance  cone.  The  50-  and  30-Q  designs  touch  the 
enclosing  sphere  at  the  top.  The  minimum  (^-factor  antenna,  32  Q,  touches  the  enclosing  sphere  on 
the  side.  The  minimum  (^-factor  antenna,  32-Q  input  impedance,  designed  with  this  simple  model, 
does  not  fill  the  top  of  the  sphere.  A  more  complex  model  would  fill  the  top  of  the  sphere  and  reduce 
the  antenna  (J-factor. 
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In  Section  5,  computer  simulation  technology  (CST)  is  used  to  compute  the  impedance  of  the 
50-  and  30-£2  input  impedance  antennas  and  the  minimum  (Q-factor  antenna,  32  Q,  where  the  VSWR 
is  3:1  at  43.1,  40.1,  and  39.3  MHz,  respectively.  The  VSWR  drops  very  quickly  to  a  VSWR  of 
~  1.2:1;  then  the  VSWR  is  ~  1.1:1  up  to  2.5  GHz.  In  Section  6  is  the  conclusion.  The  Quasi-static 
Antenna  Design  Algorithm  generates  three  UWB  non-spherical  bulb  shapes.  The  non-spherical  bulb 
shape  performs  as  well  as  the  spherical  bulb.  A  more  complex  conical  charge  distribution  would 
provide  a  lower  (^-factor  design  by  extending  the  length  of  the  cone  and  increasing  the  size  of  the 
bulb.  Minimizing  the  Q-factor  ratio  may  increase  the  VSWR  at  higher  frequencies,  which  could  be 
caused  by  an  edge  on  the  antenna  surface.  A  different  antenna  design  near  the  minimum  (^-factor 
ratio  could  give  a  smoother  antenna  surface,  which  would  reduce  the  reflection  from  the  edge 
without  significantly  increasing  the  (^-factor  ratio  of  the  antenna. 
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2.  CONICAL  CHARGE  DISTRIBUTION  AND 
EFFECTIVE  HEIGHT  FOR  UWB  ANTENNA 


The  effective  height  was  computed  [2,  3]  for  an  asymptotic  conical  dipole  (ACD)  and  a  linear 
charge  distribution  (LCD)  [3];  both  charge  distributions  are  on  the  z-axis  of  the  dipole.  A  linear 
combination  of  these  two  charge  distributions  can  create  a  conical  feed  point  for  a  UWB  antenna; 
however,  the  bulb  shape  is  nearly  spherical  and  similar  to  the  ICC  antenna. 

A  rotationally  symmetric  conical  shaped  charge  distribution  can  model  both  the  feed  point  and 
bulb  shape  of  the  UWB  antennas.  The  surface  charge  a(l)  is  integrated  around  a  ring  at  position  l. 

The  total  charge  on  this  ring  is  defined  as  a  simple  linear  function  on  / 

?(0  =  °-a)jL  +  a((^)  =  l  p('M0d9’  (2) 

where  Ka  is  the  cone  length,  a  is  the  radius  of  the  enclosing  sphere,  l  is  the  distance  to  the  ring,  and 
p  is  the  radius  of  the  ring.  The  fractional  cone  length  is  dimensionless,  0  <  k  <  1.  The  parameter 
a  -  0.35  gives  a  very  good  approximation  of  a  50-D  input  impedance  cone  at  the  feed  point.  The 
charge-distribution-cone  half  angle  ¥  (measured  from  the  z-axis)  and  k,  the  fractional  cone  length, 
generate  a  family  of  UWB  antenna  designs.  For  the  charge-distribution-cone  half  angle  ¥  =  0, 
Equation  (2)  reduces  to  a  linear  combination  of  the  ACD  and  LCD  with  the  charge  distribution  on  the 
z-axis. 

The  effective  height  depends  only  on  the  height  of  the  charge  on  the  z-axis.  The  radius  and  height 
of  the  ring  of  charge  is  p  -l  sin  ¥  and  z  =  l  cos  ¥,  respectively. 

9(— tyimage  ~  — 9(0 monopole’ 

the  charge  as  a  function  of  z,  is 

q(z)  =  q(z/  cos  ip) /cos  ip  , 


(3) 

(4) 


which  is  plotted  in  Figure  2. 

The  surface  charge  density  on  the  antenna  (with  a  120°  slice  removed)  and  cone  is  plotted  in 
Figure  1.  The  antenna  and  cone  have  the  same  net  charge;  however,  the  charge  on  the  antenna  is 
spread  out  over  a  larger  surface  area,  which  significantly  reduces  the  surface  charge  density  on  the 
antenna.  The  image  cone  also  reduces  the  electric  field  and  surface  charge  density  on  the  bulb's  top. 
The  image  cone  increases  the  electric  field  and  surface  charge  density  near  the  feed  point.  The 
surface  charge  on  the  cone  is 

0s  (0  =  9(0/(2ttp(0)  •  (5) 

Note:  p  — >  0  and  as  a  — >■  oo.  The  net  charge  on  the  ring  is  finite.  For  display  purposes,  the  range  of 
values  plotted  is  limited  to  as  <  10  *  min  (crs).  The  charge  distribution  on  the  monopole  arm  is 
normalized  to  give 

/»Ka 

qrotai  =  q(l)dl=  1.  (6) 

Jo 
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1 


(7) 


The  effective  height  is  easily  calculated  as 

q(l)zdl  = 


hEff  =  f 

qTotal  Jo 


qTotal 


L 


Ka  cos  xb  /  /  /  \ 

q(z/  COS  xp) 


zdz. 


COS  V 


This  reduces  to 


hEff  (k.  %l>) 


x  na  2 Ka 

(1  -  a)  — Ta¬ 


cos  v- 


The  radiation  resistance  is  calculated  from  the  standard  equation 


(8) 


BRad(K,y)  =  40k2h2Eff,  (9) 

where  k  =  2jt/X  and  X  is  the  wavelength.  Figure  3  is  a  plot  of  RRad  (k,  y/)  /  (ka)2  as  a  function  of  the 
cone  geometry  y  and  k.  RRad  (k,  y/)/(ka)2  is  independent  of  frequency  and  the  enclosing  sphere's 
radius  a. 


Net  Charge/(m)  on  z-axis 


Figure  2.  The  50-Q  input  impedance  monopole  and  image  charge  distribution  on  a  ring  is  integrated 
to  give  qz  net  charge/m  at  each  point  on  the  z-axis.  The  feed  point  is  at  z  =  0. 
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Charge  distribution  half  cone  angle  -  y  (rad) 


Resistance/(ka)2 


Fractional  charge  distribution  cone  length  -  k 


Figure  3.  FtRad  (x,  y)/(kaf  as  a  function  of  the  charge  distribution  geometry  k  and  y, 
which  is  computed  from  the  effective  height  of  the  charge  distribution. 
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3.  CAPACITANCE  FOR  UWB  ANTENNA  DESIGN 


The  capacitance  is  computed  from  the  net  charge  on  the  dipole  arm  and  the  maximum  scalar 
potential  on  the  enclosing  sphere.  The  scalar  potential  for  the  cone  charge  distribution  is  rotationally 
symmetric  (z,  p)  where  p  is  the  radial  distance  from  the  2-axis  and  2  is  the  position  on  the  2-axis. 
The  rotational  symmetry  reduces  the  scalar  potential  calculation  to  the  x-z  plane. 

The  cone-charge-distribution  scalar  potential  is  computed  by  replacing  the  cone's  surface  with  2N 
wires  running  from  the  charge  distribution  bi-cone's  bottom  edge  to  the  charge  distribution  bi-cone's 
top  edge.  Figure  4  shows  how  wires  are  used  to  approximate  the  conical  charge  distribution,  N  =  2 
and  N  =  6.  The  blue  solid  and  dash  lines  represent  the  N  =  2  wire  approximation.  The  blue  solid  and 
dash  lines  are  symmetric  about  the  x-z  plane.  The  N  =  6  wire  approximation  uses  the  N  =  2  results. 
The  red  dash-dot  and  dot  lines  represent  the  eight  additional  wires  required  for  theN  -  6  wire 
approximation.  The  red  dash-dot  and  dot  lines  are  also  symmetric  about  the  x-z  plane.  Only  half  of 
the  wires  (solid  and  dash-dot)  are  needed  in  the  numerical  method.  The  cylindrical  version  of  this 
numerical  approximation  is  discussed  in  detail  in  [16];  the  sequence  of  wire  approximations  N,  3N, 
9/V,...,3"W  etc.  reuses  all  of  the  previous  calculations,  which  was  suggested  by  Lance  Koyoma. 


The  N-2  and  N=6  Wire  Approximation 
of  the  Potential 


0.8-1 

0.6- 

0.4- 

0.2- 

0- 

-0.2- 

-0.4- 

-0.6- 

-0.8  - 
-0.8 


wire  1 
wire  2 
wire  1 
wire  4 
wire  3 
wire  6 
wire  2S 
wire  IS 
wire  6S 
wire  3$ 
wire  4S 
wire  IS 


”i - 1 - 1 - r - 1 - r 

-0.6  -0.4  -0.2  0  0.2  04 

y  axis 


Figure  4.  The  2-wire  (blue)  and  6-wire  (red)  numerical  approximation 
of  the  cone  charge  distribution. 
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The  scalar  potential  for  the  ACD  was  previously  computed  [2,  3]  (a  constant  line  and  image  charge 
on  the  z-axis).  The  following  equation  includes  the  image  monopole: 


<1  ’ACD(z.p) 


1 


47 TtfyKCl 


In 


(1  "T  )  (1  —  &i) 
.(1  —  3m)  (1  +  &i) . 


(10) 


where  8m  —  Ka/{Rf  +  Rt)  is  for  the  monopole  and  8t  —  Ka/{Rf  +  Rb)  is  for  the  image  monopole 
with  a  unit  net  charge  on  each  arm.  In  this  case,  the  quantities  Rt,  Rf,  and  Rb  are  the  distance  from  the 
field  point  (z,  p)  to  the  top  of  the  wire,  feed  point,  and  the  bottom  of  the  wire,  respectively.  The  scalar 
potential  for  the  linear  charge  distribution  (LCD)  on  the  z-axis,  2P1(z'/Ka)/(Ka),  was  previously 
calculated  [3]: 


/na 

— 

Ka  Ka 


2 Pi  (*7«a) 


,dz'  = 


9- 


+  P2 


(«a) 


In 


1  +T 

1  —  r 


(ii) 


where  x  -  2ica/(Rt  +  Rb)  and  z'  is  z-axis  integration  variable  for  the  charge  distribution  and  the  net 
charge  on  the  arm  is 

9  rKa 

1  =  —  /  Pi  C z'/na )  di.  (12) 

The  next  step  is  to  compute  R„  Rf,  and  Rb  for  the  wire  on  the  cone.  The  wire  is  rotated  P  radians 
around  the  rc-axis;  P  -  0  on  the  x-z  plane.  The  cone  half  angle  is  W  radians  from  the  z-axis.  The 
wire  end  point  at  the  top  cone  edge  is 


Z-t  = 

aK  cos  xjj 

(13) 

yt  = 

aK  sin  xp  sin 

(14) 

xt  = 

aK  sin  xpcos  <p. 

(15) 

The  same  wire  end  point  on  the  bottom  cone  edge  is 


Zb  - 

— aK  cos  xp 

(16) 

lib  — 

—aK  sin  xp  sin  <p 

(17) 

xb  = 

—aK  sin  xp  cos  <p. 

(18) 

The  distance  from  the  top  cone  edge  to  the  point  (z,  p)  (in  the  x-z  plane)  is 


Rt  (<p)  =  \j [s  -  «a  cos  V']2  +  [p  - 


aK  sm  vj  cos  ^p\ 


+  [< 


an  sin  v  sin 


(19) 


The  distance  from  p,  z  (in  the  x-z  plane)  to  the  feed  point  of  the  monopole  is 

Rf  =  \/s2  +  /o2  • 


(20) 


The  distance  from  p,  z  (in  the  x-z  plane)  to  the  bottom  cone  edge  (image  monopole)  is 

Rb  {<?)  =  yj\z  +  Ka  cos  +  [p  +  an  sin  xp  cos  y?]2  4-  [aK  sin  xp  sin  <p]2 


(21) 
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The  variable  8m(<p)  =  a/c/[/?t(<p)  +  Rf]  is  for  the  monopole  and  8i(cp')  =  a.K/[Rb((p)  +  fy] 
is  for  the  image  monopole.  For  the  N  wire  case, 


‘Pj  =  (j  —  1/2)  n/N 


with  1  <  j  <  N . 


The  scalar  potential  is 

[l  +  <5m  (y>,)]  [l  -  6j  (<?,-)] 

[1-^m  fo)]  [1  +  8i  fo)] 


(22) 


(23) 


The  wires  with  fj  =  -  (j  ~i/2)-/v  give  the  exact  same  equation  as  above  and  are  not  included  in 
the  calculation. 


The  LCD  scalar  potential,  Equation  (11),  assumes  the  line  charge  is  on  the  z-axis  with  a  unit 
vector  Zj  =  1,  Xj  =  0,  y;-  =  0  .  The  dot  product  of  the  unit  vector  and  the  field  position  is  z.  For  the 
cone  case,  the  unit  vector  for  the  wire  in  direction  cpj  and  ip  is  z;-  =  cosi|/,  Xj  —  sin(xp)cos(pj,  yj  — 
sin(xp)sin(pj. 

The  field  position  is  x  —  p  and  y  —  0.  The  projection  Zj  of  the  field  position  on  the  cpj  wire  is 
calculated  with  the  dot  product: 


Zj  =  z  cos  xb  +  p  sin  xp  cos  <p3 . 


(24) 


The  scalar  potential  is 


<I> 


LCD 

N 


(*>/>) 


7 = N 

2  Zi 


4t tcq  (h 


Ef 


i=i 


(25) 


In  the  general  case,  the  error  in  the  N  wire  numerical  approximation  is 


*£rror  (*,/>)  =  |4>3at  (*,/>)-  <I>v  (*,  p)  |-  (26) 

The  ‘N.v  U, p)includes  the  wires  from  ‘Lv  {?-,  p)  plus  two  more  wires  on  either  side  of  the  N  wires.  The 
numerical  method  error  fffrror  (z,p)  is  calculated  to  give  ^W  or  (z,  p)  /l^s.v  (z,  p)  \  <  to-7.  The  exact 
error  ^fjvror  (Z,P>  can  only  be  evaluated  after  T9  V  (z,  p )  is  computed;  however,  the  error 

$3N°r  (*>  p)  «  4>frror  (»,  p)  I  (27) 


is  much  smaller.  The  leading  error  term  can  be  deduced  by  symmetry  and  spherical  harmonics. 

The  scalar  potential  and  cone  are  both  rotationally  symmetric;  however,  the  error  in  the  solution 
has  the  same  rotational  symmetry  of  the  2N  wire  numerical  approximation.  The  error  is  found  with 
the  spherical  harmonic  expansion  of  rP).  The  wires  are  equally  spaced  with  &4>  =  ,TAY  steps 

in  rotational  angle  and  they  are  odd  in  z. 

The  spherical  harmonic  term,  Tm  =  \/2l  + 1  (/  -  ~  mi'-PR  (cos  0)eim  contributes  to  the  expansion 

only  if  it  has  the  same  rotational  symmetry  as  the  wires  m  =  0.  2 N,  AN,  etc.  Note: 
mA<P  =  2N  * 71 /N  =  2w.  The  associated  Legendre  polynomials  must  be  odd  pr(cos  9)  =  -i7*(cos(ir-  9)\ 
this  limits  /  to  odd  values.  The  solution  is  the  rotationally  symmetric  part  with  no  e,m<>  dependence; 
the  only  spherical  harmonic  expansions  contributing  to  the  solution  are  f2P+i  o/r2”+-  where  0  <  p  <  oo. 
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The  first  error  term  in  the  spherical  harmonic  expansion  is  T2,v+i,2,v/r2A_2.  This  term  is  rotationally 
symmetric  in  ttN  steps  and  odd  in  z.  The  first  error  term  in  the  spherical  harmonic  expansion  of 
&3N(z,p)  is  Y6N+16N/r6N+2.  At  large  distances,  the  largest  error  is  proportional  to 

$frrror  (z,p)  oc  \Y2N+h2N/r2N+2\.  (28) 

This  numerical  method  converges  very  fast  with  the  sequence  N,  3 N,  9 N,  etc.  The  ACD  and  LCD 
scalar  potentials  are  combined  with  the  coefficients  from  Equation  (2): 

*/TiW)  =  (!-<*)  <CD{*,P)  +  ^lncd(z,p).  (29) 

The  capacitance  is  computed  from 

TmLB  («»  =  max  B  (  ^Sphere  ■  p Sphere )]  ’  (30) 

where  zsPhere  =  a  cos  9.  P  sphere  =  asm  9  and  9  is  measured  from  the  z-axis  and  0  <  6  <  n/2. 

The  capacitance  is 


c"*fc«=WC*(>.«'  <31> 

Note:  the  charge  distribution  is  defined  to  give  tfTotai  =  1,  Equation  (6).  Figure  5  is  the 
capacitance/a  as  a  function  of  the  cone  geometry  ip  and  k.  The  capacitance  is  scaled  by  1/a  to 
make  the  design  curve  independent  of  the  enclosing  sphere’s  radius  a.  For  a  fixed  ip,  increasing 
k  moves  the  charge  distribution  cone  closer  to  the  enclosing  sphere;  this  increases  the  peak 
scalar  potential  and  decreases  capacitance.  For  a  fixed  k,  the  larger  half  cone  angle  i//  increases 
the  distance  between  the  enclosing  sphere  and  the  charge  distribution  on  the  cone.  This 
decreases  the  peak  scalar  potential  on  the  enclosing  sphere.  The  negative  charge  on  the  image 
cone  is  also  closer  to  the  top  part  of  the  enclosing  sphere;  this  also  reduces  the  peak  scalar 
potential.  The  peak  scalar  potential  on  the  enclosing  sphere  is  significantly  reduced  for  the 
higher  half  cone  angles;  this  increases  the  capacitance.  In  Figure  5,  the  distance  between 
contours  is  smaller  for  the  larger  values  of  ip. 
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Charge  distribution  half  cone  angle  -  y  (rad) 


0.8 


Capacitance/a  (picoF/m) 


Figure  5.  Capacitance  design  curve  as  function  of  the  charge  distribution  geometry  k  and  ijj,  which  is 
computed  from  the  potential  on  the  antenna  and  net  charge  on  dipole  arm. 
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4.  UWB  ANTENNA  Q-FACTOR  RATIO 

The  capacitance  and  radiation  resistance  are  combined  to  compute  the  Qfactor  ratio  as  a  function 
of  fractional  charge-distribution-cone  length  and  charge-distribution-cone  half  angle  y.  Only  a  subset 
of  solutions  representing  an  antenna  with  a  conical  feed  point  will  be  considered.  Other  solutions 
with  non-conical  feed  points  may  not  be  UWB.  The  value  a  =  0.35  is  selected  for  the  50-  Q  input 
impedance  antenna  design;  it  is  a  good  approximation  of  the  feed  point  cone  for  the  other  design  in 
the  paper.  The  accuracy  of  the  feed  point  cone  is  critical  to  the  calculation  of  the  capacitance/a, 
radiation  resistance/i/a)2,  and  Qfactor  ratio.  In  the  following  plots,  the  antenna  design  curves  are 
independent  of  frequency  and  enclosing  sphere's  radius  a.  The  Qfactor  ratio  is  valid  in  the  limit  as 
frequencies/^  0. 

Figure  6  plots  Qfactor  ratio,  Q  *(kaf,  as  a  function  of  the  cone  geometry  k  and  t/A  For  small  <//, 
the  minimum  (/-factor  ratio  has  moderate  ><  values.  For  large  VO  the  (/-factor  ratio  has  high  k  values. 
The  transition  region  is  expanded  in  Figure  7;  there  is  an  absolute  minimum  and  a  local  minimum  for 
some  values  of »//.  The  absolute  and  a  local  minimum  are  the  same  at  »//  =  0.19481 ;  the  absolute 
minimum  shifts  from  k  ~  0.84  to  n  ~  0.98. 


Q-factor  Ratio  Q*<ka)3 


Figure  6.  The  Q-factor  ratio  design  curve  as  function  of  the  charge  distribution  geometry  k  and  y. 
The  Q-factor  ratio  design  curve  is  computed  from  the  radiation  resistance  and  capacitance. 
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Q-factor  Ratio  Q*(ka)3 


Figure  7.  The  Q-factor  radio  has  an  absolute  minimum  and  a  local  minimum  for  a  range  of  cone 
geometries  k  and  y.  As  y  increases,  the  absolute  minimum  shifts  from  k  ~0.84  to  k  ~0.98. 

Figure  8  plots  the  minimum  Q-factor  ratio,  Q  *  ( ka )3,  as  a  function  of  ¥  in  radians;  the 
capacitance,  C/a;  radiation  resistance,  RRaJ(ka)2;  and  10k  for  this  solution  are  also  plotted.  At  the 
discontinuity,  Q-factor  ratio  has  the  same  values  at  two  point  with  different  values  of  k.  Below  the 
discontinuity,  the  increasing  k  in  turn  increases  RRa</(ka)2  and  also  increases  the  scalar  potential  on 
the  sphere,  which  reduces  the  capacitance  of  the  antenna.  Above  the  discontinuity,  k  is  almost  a 
constant,  RRa/(kaf  decreases  with  the  cos2¥.  The  capacitance  increases  with  the  charge-distribution- 
cone  half  angle  ¥.  The  surface  charge  density  is  spread  over  a  much  larger  surface  area,  which 
reduces  the  peak  scalar  potential  on  the  enclosing  sphere  and  increases  the  capacitance.  For  ¥  > 

0.177,  the  minimum  Qfactor  antenna  no  longer  touches  the  top  of  the  enclosing  sphere.  For  charge- 
cone  half  angles  0.195  <¥<  0.715  range,  k  is  constant. 

The  first  row  in  Table  1  is  the  minimum  Qfactor  ratio  for  the  line  charge  distribution.  The  last  two 
rows  in  Table  1  have  similar  input  impedance  but  very  different  shapes.  Figure  9  plots  the  line- 
charge-distribution  Qfactor  ratio  design  curve  as  a  function  of  k.  Figure  10  plots  the  antenna  shape 
for  the  ICC  antenna  and  the  minimum  Qfactor-ratio,  line-charge-distribution  design.  The  line- 
charge-distribution  bulb  shape  is  almost  the  same  as  the  ICC  antenna,  but  the  ICC  feed  point  cone  is 
much  smaller.  The  50-  Q  input  impedance  ICC  antenna  will  be  used  as  the  best  existing  solution. 
Simpson,  Pavlovic,  and  Olcan  [13]  did  not  evaluate  the  ICC  antenna  for  different  cone  half  angles; 
the  50-  Q  input  impedance  ICC  antenna  may  not  be  the  minimum  Q  -factor  solution. 
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Minimum  Q  Factor  Ratio  Design  Parameters  as  a 
Function  of  Charge  Distribution  Cone  Angle 


Charge  distribution  cone 
half  angle  -  y  (rad) 


Figure  8.  The  minimum  Q-factor  ratio,  capacitance/a,  10  k,  and  resistance/^)2  as  a  function  of  <//. 

At  the  discontinuity,  the  Q-factor  ratio  has  the  same  value  at  two  different  points;  see  Figure  7  and 
Figure  1 1 . 

The  second  row  in  Table  1  is  the  maximum  resistance  design.  Figure  1 1  plots  Qfactor-ratio  design 
curves  for  i//  =  0.1946,  y/  =  0.19481,  and  y/  =  0.195.  The  minimum  Q  shifts  from  k  ~  0.84  on  the  /// 

=  0.1946  design  curve  to  k  ~  0.98  on  the  i//  =  0.195  design  curve.  The  design  curve  for  y/  =  0.19481 
has  two  equal  minimum  (^-factor-ratios.  The  maximum  Rruci  design  does  not  give  a  significant 
decrease  in  Q-factor  ratio  and  will  not  be  numerically  modeled. 

Row  3  in  Table  1  is  50-  T2  input  impedance  design.  The  parameters  are  selected  to  touch  the 
enclosing  sphere  at  the  top.  The  50- T2  input  impedance  design  curve.  Figure  12,  shows  the  minimum 
Qfactor  ratio  and  the  position  of  the  design.  Figure  13  is  the  50- Q  input  impedance  design  cross 
section  that  fills  the  top  of  the  sphere.  Figure  14  shows  that  the  antenna  design  is  a  very  good 
approximation  of  the  50-  Q  input  impedance  cone. 
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The  last  two  rows  in  Table  1  have  similar  input  impedance  but  very  different  shapes.  A  30-input 
impedance  design  touches  the  top  of  the  enclosing  sphere.  On  the  other  hand,  the  32-minimum 
Q-factor-ratio  design  touches  the  side  of  the  enclosing  sphere.  Figure  15  shows  the  two  design  curves 
for  fixed  \\i  =  0.5296  and  \|/  =  0.653.  There  is  a  significant  difference  on  the  Q-factor  ratios  of  the  two 
antennas.  Figure  16  shows  there  is  a  significant  difference  in  the  antenna  shape.  The  next  section 
shows  that  the  32-input  impedance  design  reduces  the  lowest  operational  frequency  without 
increasing  the  VSWR  at  higher  frequencies. 

Table  1.  Antenna  design  parameters:  radiation  resistance/^/,  capacitance,  and  (^-factor  ratio. 


Parameter 

¥ 

^Rac/fe)2 

da 

Q*(ka)3 

Line  charge 

0 

0.7150 

6.3747 

96.6077  pF 

5.413 

Max  RRad 

0.1948 

0.9789 

11.5 

57.45  pF 

5.045 

50  Q 

0.3408 

0.940 

9.7871 

84.9754  pF 

4.008 

30  Q 

0.5296 

0.8711 

7.047 

1 27.805  pF 

3.701 

Min.  Q  factor  ratio  32  Q 

0.653 

0.9809 

7.5688 

134.890  pF 

3.265 

Q-factor-Ratio  for 
Line  Charge  Distribution 


length  -  k 

Figure  9.  (^factor  ratio  for  the  z-axis  line  charge  distribution  as  a  function  of  k  with  ¥  =  0. 


16 


Minimum  Q  Line  Charge  Distribution  and  ICC 


Figure  10.  The  50-0  input  impedance  ICC  and  the  minimum  Q line  charge 
distribution  designs  (V  =  0). 


Q-factor-Ratio  for  Max  R„  . 

Rad. 


Figure  1 1 .  The  k  of  the  minimum  (^shifts  from  0.84  to  0.979.  This  significantly 
increases  the  FtRad  without  changing  Q. 
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Q-factor-Ratio  for  the  50  Q 


Fractional  charge  distribution  cone  length  -  k 


Figure  12.  The  (^-factor  ratio  design  curve  for  charge  distribution  cone 
geometry  w  =  0.3408  and  k.  The  design,  marked  with  touches  the 
top  of  the  sphere. 


A  1 -meter  high  50  Q  Antenna  Design 


Figure  13.  The  enclosing  sphere,  antenna  shape,  feed  point  cone  angle, 
and  charge  distribution  cone. 
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The  z-axis  height  above  the  feed  point  (m) 


50Q  Feed  Point 


near  the  feed  point  (m) 


Figure  14.  The  antenna  shape  is  a  very  good  approximation  of  the  50- Q  cone. 
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Minimum  Q-Factor-Ratio 


Fractional  charge  distribution  cone  length  -  k 

Figure  1 5.  The  Q-factor  ratio  is  significantly  higher  for  the  antenna  that  touched  the  top  of  the 
enclosing  sphere,  V  =  0.5296.  The  design  that  touched  the  side  of  the  enclosing  sphere  has  a  lower 
Q-factor  ratio,  V  =  0.653. 


Two  Antenna  with  Similar  Q-factor  Ratios 


Figure  16.  The  green  antenna  has  a  lower  Q-factor  ratio  than  the  blue  antenna.  The  higher 
curvature  on  the  green  antenna  does  not  influence  the  VSWR  of  the  antenna,  see  Figure 
27. 
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5.  NUMERICAL  RESULTS 


The  numerical  results  are  calculated  with  the  T-solver  in  Computer  Simulation  Technology  (CST) 
Microwave  Studio®.  The  numerical  model  uses  a  coaxial  cable  (source)  embedded  in  a  thin  cone 
6  mm  high  with  a  diameter  of  400  mm  on  the  ground  plane.  For  the  50-  Q  coaxial  source,  the  inner 
conductor  is  6  mm  radius  wire  and  the  outer  conductor  is  13.8  mm  in  radius.  Air  is  the  insulator 
between  the  conductors.  The  coaxial  inner  conductor  extends  into  the  body  of  the  antenna  replacing 
the  point  of  the  cone  with  a  cylinder.  The  problem  size  is  reduced  by  using  B  =  0  on  the  x-z  and  y-z 
symmetry  plane.  Energy  adaptation  was  used  to  obtain  a  series  of  mesh  refinements.  In  the  initial 
work,  the  antennas  were  selected  to  touch  the  top  of  the  sphere,  p  -  0  and  z  -  1. 

Figure  17  is  the  cross  section  of  the  CST  antenna  model  for  the  50- Q  input  impedance  antenna. 
Figure  18  is  the  Smith  chart  and  Figure  19  is  the  VSWR  plot  for  the  frequency  range  30  MHz  to 
2.5  GHz.  The  VSWR  is  2:1  at  54.813  MHz;  this  is  an  8.6%  reduction  in  size  compared  to  the  results 
of  Simpson,  Pavlovic,  and  Olcan  [13],  The  VSWR  is  less  than  1.2:1  for  nearly  all  of  the  frequency 
range.  This  is  an  unexpected  result. 


Figure  17.  A  cross  section  of  the  CST  50- Q  input  impedance  1 -meter-high  monopole. 


Figure  18.  The  Smith  chart  for  the  1 -meter-high  50- Q  input 
impedance  monopole. 
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Figure  1 9.  The  VSWR  for  the  1  -meter-high  50-  Q  input  impedance  monopole. 


The  Q-factor  and  lowest  operating  frequency  can  be  reduced  by  increasing  charge-distribution- 
cone  half  angle  V;  this  lowers  the  impedance  of  the  feed-point  cone  and  increases  the  capacitance. 
The  numerical  model  uses  a  30-  Q  coaxial  source  embedded  in  a  thin-cone  3  mm  high  with  a 
diameter  of  400  mm  on  the  ground  plane.  For  the  30-  Q  coaxial  source,  the  inner  conductor  is  3  mm 
in  radius  wire  and  the  outer  conductor  is  5.001  mm  in  radius.  Air  is  the  insulator  between  the 
conductors.  The  coaxial  inner  conductor  extends  into  the  body  of  the  antenna  replacing  the  point  of 
the  cone  with  a  cylinder.  Energy  adaptation  was  used  to  obtain  a  series  of  mesh  sizes.  This  antenna 
touches  the  top  of  the  sphere,  p  =  0  and  z  =  I . 

Figure  20  is  the  cross  section  of  the  CST  antenna  model  for  the  30-  Q  input  impedance  antenna. 
Figure  21  is  the  Smith  chart  and  Figure  22  is  the  magnitude  of  the  reflection  coefficient,  T  for  the 
frequency  range  20  MHz  to  2.5  GHz.  The  VSWR  is  2: 1  at  48.9  MHz  or  X/6.134.  This  is  almost  a 
18.5%  reduction  in  size  compared  to  Simpson,  Pavlovic,  and  Olcan’s  [13]  results.  The  reflection 
coefficient  ranges  from  1/10  to  1/20  for  nearly  all  of  the  frequency  range  (VSWR  1.2  to  1.1).  Figure 
23  is  the  CST  source  pulse  and  Figure  24  is  the  CST  reflected  pulse.  Note:  the  scale  on  the  reflected 
pulse  is  4%  of  source  pulse.  At  the  feed  point,  the  reflection  is  3.9%  of  the  incident  pulse. 
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Figure  20.  Thel -meter-high  30- Q  input  impedance  UWB  antenna  design  that  touches 
the  top  of  the  sphere. 


Figure  21 .  The  Smith  chart  for  30- Q  input  impedance  antenna  that  touches  the  top  of  the 
sphere. 
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S-Para meter  [Magnitude] 


Frequency  /  MHz 


Figure  22.  The  CST  frequency  domain  reflection  coefficient  for  the  30- Q  input  impedance  antenna 
that  touches  the  top  of  the  sphere.  The  VSWR  4:1 , 3:1 ,  and  2:1  are  labeled  as  1 , 2,  and  3, 
respectively. 


Figure  23.  The  CST  source  pulse  for  30- Q  input  impedance  antenna  that  touches  the  top  of  the 
sphere. 
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Figure  24.  The  CST  pulse  reflected  from  the  feed  point  for  30- Q  input  impedance  antenna  that 
touches  the  top  of  the  sphere. 


The  final  design  places  no  restriction  on  the  antenna  shape.  The  top  of  the  antenna  no  longer 
touches  the  top  of  the  sphere.  Figure  25  is  the  cross  section  of  the  CST  antenna  model  for  the  32- Q 
input  impedance  antenna.  Figure  26  is  the  Smith  chart  and  Figure  27  is  the  VSWR  plot  for  the 
frequency  range  of  20  MHz  to  2.5  GHz.  The  VSWR  is  2: 1  at  48.9  MHz,  which  is  slightly  larger  than 
the  previous  result.  The  VSWR  is  less  than  1.2:1  for  nearly  all  of  the  frequency  range.  Figures 
28  and  29  are  the  CST  source  and  reflected  pulse.  Note:  the  scale  on  the  reflected  pulse  is  4%  of 
source  pulse.  At  the  feed  point,  the  reflection  is  3.9%  of  the  incident  pulse.  In  Figure  16,  the 
minimum  Q antenna  design  has  a  rapid  change  in  curvature  near  the  enclosing  sphere.  One  would 
expect  RF  to  be  reflected  from  this  region,  but  the  RF  is  not  reflected,  which  implies  that  other 
changes  in  shape  can  be  made  without  increasing  the  VSWR.  The  antennas  capacitance  and  radiation 
resistance  can  both  be  increased  and  Q-factor  decreased  by  filling  the  top  of  the  sphere.  A  more 
complex  charge  distribution  model  can  be  used  to  extend  the  cone  length,  fill  the  top  of  the  sphere, 
and  significantly  reduce  the  Q-factor. 

Table  2  summarizes  the  VSWR  numbers  for  the  three  designs.  The  frequencies  for  3  to  1  and  4  to 
1  VSWR  are  steadily  decreasing  with  lower  Q.  Hujanen,  Holmberg,  and  Sten  [17]  calculated  the 
smallest  possible  ka  value  for  the  UWB  antenna  size  and  VSWR.  This  analytic  solution  assumes  a 
dipole  mode  with  an  infinite  lossless  matching  network.  The  UWB  antenna  size  is  ka  0.46833  or  a 
=  Amax/13.4  for  a  3: 1  VSWR.  They  pointed  out  that  an  antenna  can  be  “pre-tuned”  and  not  require  a 
matching  network.  Their  derivation  does  not  indicate  the  size  of  a  pre-tuned  antenna.  Yang,  Davis, 
and  Stutzman  [18]  designed  an  antenna  with  a  10:1  bandwidth  and  an  antenna  size  a  =  Xmax/8. 
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Table  2.  Lowest  operating  Frequency  for  Different  VSWRs. 


Design 

Q  *  ( ka )3 

p  =  1/3 

2:1  VSWR 

3:1  VSWR 

4:1  VSWR 

50  Q 

4.008 

1.149  ka 

54.81  MHz 

47.18  MHz 

43.08  MHz 

30  Q 

3.701 

1 .006  ka 

48.005  MHz 

43.23  MHz 

40.05  MHz 

Min.  Ofactor  ratio  32  Q 

3.265 

1 .025  ka 

48.90  MHz 

42.67  MHz 

39.25  MHz 

Figure  25.  Minimum  Q-factor  design  computed  from  the  QSADA. 
The  antenna  does  not  fill  the  top  of  the  sphere. 


Figure  26.  Smith  chart  for  the  minimum  Q-factor  design  computed  from  the  QSADA.  The  VSWR 
4:1, 3:1, 2:1,  and  1.14:1  are  labeled  by  1, 2,  3,  and  4,  respectively. 
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Figure  27.  VSWR  for  the  minimum  ^factor  design  CST  model. 
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Figure  28.  The  CST  source  pulse  (red)  and  reflection  (green)  for  the  minimum  Q-factor  design 
computed  from  the  QSADA. 
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Figure  29.  The  CST  reflection  from  the  feed  point  for  the  minimum  0-factor  design  computed  from 
the  QSADA. 
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6.  CONCLUSION 


The  Quasi-static  Antenna  Design  Algorithm  is  a  general  method  for  designing  cylindrically 
symmetric  antenna.  This  report  covers  the  UWB  antenna  design  and  [6]  discusses  the  electrically 
small  antenna  design.  The  general  top  load  [6]  can  be  added  to  the  conical  charge  distribution  to 
model  more  complex  charge  distributions. 

The  Quasi-static  Antenna  Design  Algorithm  is  a  physics-based  approach  for  UWB  antenna 
designs.  The  cone  with  a  bulb  shape  is  used  as  a  UWB  antenna  model.  Minimizing  the  Q-factor  ratio 
for  this  model  yields  cone  and  bulb  shapes  with  improved  low-frequency  performance.  The  lowest 
(Q-factor  ratio  is  a  32-  Q  input  impedance  UWB  antenna  that  fills  most  of  the  enclosing  sphere's 
volume.  The  30-  Q  input  impedance  and  minimum  Q-factor  ratio  designs  have  very  different  shapes 
with  similar  (Q-factor  ratios  and  very  low  VSWR,  as  shown  in  Figure  16.  The  shape  difference  should 
reflect  RF  back  to  the  feed  point;  this  is  an  unexpected  result.  One  would  expect  that  the  VSWR  to 
increase  as  the  operating  frequency  is  reduced.  The  lowest  operating  frequency  of  the  UWB  antenna 
should  give  a  uniform  VSWR  over  the  entire  band.  Any  further  reduction  in  operating  frequency 
would  reduce  the  antenna  bandwidth. 

The  above  model  uses  a  linear  charge  distribution  on  the  cone.  The  quasi-static  antenna  design 
algorithm  can  use  a  sequence  of  more  complex  charge  distribution  models:  quadratic,  cubic,  etc.,  on 
the  cone.  Each  new  charge  distribution  should  extend  the  length  on  the  cone  and  the  size  of  the  bulb, 
which  moves  the  charge  distribution  to  a  higher  position  on  the  antenna,  increasing  the  effective 
height  and  radiation  resistance.  The  larger  cone  and  bulb  will  also  reduce  the  surface-charge  density 
on  the  antenna,  which  reduces  the  electric  field  around  the  antenna  and  increases  the  capacitance. 
Both  factors  reduced  the  Q.  Each  minimum  (Q-factor  design  requires  numerical  modeling  to  evaluate 
the  impedance  and  the  VSWR.  The  quasi-static  antenna  design  algorithm  will  reduce  the  lowest 
operating  frequency  of  the  antenna;  however,  the  new  antenna  shape  could  adversely  impact  the  high 
frequency  response  of  the  antenna.  The  Q-factor  design  curves  allow  the  antenna  shape  to  be 
modified  by  selecting  a  new  design  near  the  minimum  Q factor  ratio.  This  could  round  the  bulb 
shape  and  reduce  VSWR.  The  Quasi-static  Antenna  Design  Algorithm  should  be  able  to  design 
significantly  smaller  UWB  antennas. 
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